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The vertical radiation force on an absorbing micrometer-sized dielectric sphere situated in an 
evanescent field is calculated, using electromagnetic wave theory. The present work is a continuation 
of an earlier paper [E. Almaas and I. Brevik, J. Opt. Soc. Am. B 12, 2429 (1995)], in which both the 
horizontal and the vertical radiation forces were calculated with the constraint that the sphere was 
non-absorbing. Whereas the horizontal force can be well accounted for within this constraint, there 
is no possibility to describe the repulsiveness of the vertical force, so distinctly demonstrated in the 
Kawata-Sugiura experiment [Opt. Lett. 17, 772 (1992)], unless a departure from the theory of pure 
non- dispersive dielectrics is made in some way. Introduction of absorption, i. e. a complex refractive 
index, is one natural way of generalizing the previous theory. We work out general expressions for 
the vertical force for this case and illustrate the calculations by numerical computations. It turns out 
that, when applied to the Kawata-Sugiura case, the repulsive radiation force caused by absorption 
is not strong enough to account for the actual lifting of the polystyrene latex or glass spheres. 
The physical reason for the experimental outcome is, in this case, most probably the presence of 
surfactants making the surface of the spheres partially conducting. 

PACS numbers: 03.50.De, 42.25.-p, 42.50.Vk, 42.55.-f 



I. INTRODUCTION 

The outcome of the Kawata-Sugiura levitation exper- 
iment from 1992 Q] is in some ways still surprising. The 
authors examined the radiation force on a micrometer- 
sized spherical dielectric particle in the evanescent field 
produced by a laser beam of moderate power, P — 150 
mW. As one would expect, the particle moved along the 
horizontal plate at a speed of a few micrometers per sec- 
ond, as a result of the tunneled photons in the evanescent 
field. The surprising point is, however: why does the par- 
ticle lift from the surface? As reported in Ref. "the 
particle is forced to float from the substrate surface and 
to slide along the surface." For a non-magnetic and non- 
absorbing dielectric particle, it is well known that the 
electromagnetic volume force density f is 



1 



(1) 



e being the permittivity such that f = in the homoge- 
neous and isotropic interior, and different from zero only 
in the particle's boundary layer. The evanescent field is 
known to vary as exp[— /?(a; -I- h)] in the vertical direction 
(see Eq. ()14|l below) , and is significant only very close to 
the substrate. In practice, this field interacts with those 
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parts of the particle that are situated closest to the sub- 
strate, and it gives rise to a surface force that, according 
to Eq. JQ) , has to act in the downward direction (assum- 
ing that e > eo). It is simply impossible to account for 
a lift force on the particle, using Eq. This discrep- 
ancy between electromagnetic theory for dielectrics and 
the Kawata-Sugiura experiment is therefore, as men- 
tioned above, somewhat surprising at first sight, since 
Eq. Q has generally been proven to be invaluable in a 
multitude of cases in electromagnetism and optics. 

Consequently, we have to conclude that the explana- 
tion for the observed lift of the particle has to lie in the 
presence of an additional force, different from the one 
given by Eq. . The following two possibilities come to 
mind: 

(1) The particle has a complex permittivity e, and is 
therefore exposed to a repulsive (or positive lift) force 
caused by the absorption of radiation in the interior. 

(2) Another natural possibility is that surface effects 
may have played a significant role in the experiments. 
The presence of impurities on the surface, or the pres- 
ence of a liquid of a film of adsorbed material, may have 
a considerable influence on the electrical conductivity of 
the surface region, and can thus enhance multiple scat- 
tering effects. This allows for an analysis of the problem 
in terms of a multiple scattering formalism. Adopting 
such a picture, the vertical force on the particle neces- 
sarily has to be repulsive; the photons are bouncing off 
the spherical surface and are, accordingly, transferring an 
impulse in the upward direction. 

In the present paper, we will be studying option (1) 
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above, thus allowing for a complex e, which in turn im- 
plies a non- vanishing effective conductivity a in the inte- 
rior of the particle. As far as we know, such a study has 
not previously been undertaken. This study is a continu- 
ation of the work of Almaas and Brevik j^J , dealing with 
the case of a pure dielectric (real e). Whereas the analysis 
in gave satisfactory results for the horizontal radiation 
force observed in the Kawata-Sugiura experiment, it was 
impossible, as mentioned above, on the basis of Eq. 
to account for an upward directed vertical force. From 
now on, we will assume that the particle is a sphere of 
radius a, and regard the effective conductivity a as an 
input parameter. 

Before embarking on the mathematical formalism, let 
us however make some further remarks related to option 
(2) above. There is a variety of surface effects that one 
can envisage: multiple reflections, as mentioned; more- 
over capillary forces (the order of magnitude of the cap- 
illary forces may be up to 10"'' N), surface tension (be- 
cause of the capillary forces), viscous retarding forces if 
the sphere and the surface are immersed in a liquid, elec- 
trostatic forces, and finally the role of evaporation. Cf. 
in this context the survey given by Vilfan et al. 3] . Espe- 
cially the role of the size of the spherical particle seems to 
be important. If the particle is large, multiple scattering 
between the sphere and the substrate distorts the evanes- 
cent wave structure that one would have in the absence 
of the particle. Lester and Nieto-Vesperinas 3 recently 
analyzed the problem using multiple scattering methods, 
and were able to account for a positive lift force. An ad- 
vantage of this kind of method is that the proximity of the 
sphere-substrate two-body system is addressed explicitly. 
In contrast, our formalism is built on the model of scat- 
tering of an evanescent wave by an isolated sphere, and 
is therefore not expected to be very accurate for small 
separations between the sphere and the substrate. It is 
necessary to point out that the adoption of a multiple 
scattering method does not in itself imply that the lift 
force is positive, as long as the material in the sphere is 
a pure dielectric. Regardless of whether the electromag- 
netic boundary conditions at the surface of the substrate 
are taken into account or not, for a pure dielectric, the 
force on the sphere is given by Eq. and, as we have 
argued above, is necessarily downward directed. To ac- 
count for a positive lift force, we must allow for deviations 
from the simple theory of pure dielectrics, meaning that 
absorption, and consequently conductivity, are important 
effects in the theory, explicitly or implicitly. 

As already mentioned, we will not attempt to give a 
quantitative description of a possible adsorbed film on 
the surface; the physical reason for such a film may be 
quite complex. Problems of this sort belong to the vast 
field of physicochemical hydrodynamics. [A good source 
of information in this area is the book of Levich [a; a 
briefer treatment is found in Landau and Lifshitz |a[|. 
It is noticeable that the experiment of Vilfan et al. Q 
similarly measured a repulsive force. They obtained a 
force 300 pN on a dielectric sphere of radius 5 fim in the 



evanescent field of a 500 mW laser focused to a diame- 
ter of 5 /im. These spheres are relatively large, and are 
therefore likely to distort the evanescent field structure 
considerably. 

In the next section, we summarize the wave theoretical 
formalism describing the interaction between the evanes- 
cent field and the sphere. In Section 3, we show how the 
surface force F*'"''^ - i.e. the part of the total radiation 
force F which is associated with the boundary and in- 
dependent of the absorptive properties in the interior - 
can be written as a double sum over complex expansion 
coefficients. Section 4, the main section of our paper, 
contains the calculation of the absorptive part F'*'^'' of 
the radiation force. The results are illustrated in Figs. [21 
- |S| Our calculations incorporate both the two polariza- 
tion states of the incident beam (p and s polarization). 

Our main analytical result is given in Eq. (|66|l . for 
the non-dimensional vertical absorptive force component 
Q'^^. As one would expect physically, the result be- 
comes proportional to the effective conductivity a. For 
the polystyrene latex spheres, or glass spheres, used in 
the Kawata-Sugiura experiment, the values of a are low. 
Hence, it turns out that Q'^^ becomes too weak to ac- 
count for the floating of the sphere in the gravitational 
field. Now, our formalism has, of course, a much wider 
scope than the Kawata-Sugiura experiment. In optics, 
there are numerous cases where radiation forces interact 
with absorbing media. For instance, the current develop- 
ment of microelectromagnetical systems, which are used 
in a variety of applications including even measurements 
of the Casimir force, accentuates the usefulness of the 
present theory. 



II. BASIC FORMALISM 

A sketch of the geometry is shown in Fig. A 
laser beam is incident from below through a transparent 
medium, called medium 1, and hits the plane, horizontal 
surface towards the transparent medium 2 at an angle 
of incidence, 9i, which is greater than the critical angle, 
^crit = arcsin(n2/ni), for total reflection. This estab- 
lishes an evanescent field above the plane surface, which, 
in turn, is regarded as the incident field falling upon a 
compact sphere of radius a and complex refractive index 
n^. The sphere is centered at the origin with the x axis 
pointing vertically upwards. If £3 is the complex permit- 
tivity of the (non-magnetic) sphere, we will require it to 
be related to the effective conductivity, a, through the 
equation 

£3 = £3 + «cr/cj, (2) 

hence fis = i/ea/eo- 

Here, a remark is called for as regards our formal use 
of the dispersion relation. For simplicity, taking the 
medium to have only one resonance frequency, we can 



3 




h 




FIG. 1: Spherical particle of radius a, with complex refractive 
index n-j,, situated in an evanescent field with its center at a 
height h above the plane substrate. A laser beam is incident 
from below at an angle of incidence 9\ > Ocrit- Refractive 
indices in the transparent media 1 and 2 are respectively ni 
and n2. 



write the complex permittivity as 



e(w) = 63 + i 



cj(7 — iuj) ■ 



(3) 



cf., for instance, Eq. (7.56) in Jackson 0- Here the first 
term ea describes the contribution from the dipoles and 
can be taken to be constant, whereas the second term 
describes the dispersive properties of the free electrons. 
7 is the damping constant. Comparison between Eqs. (0) 
and gives for the frequency dependent term 



a{uj) 



J — iu! 



(4) 



Typically 7 ^ 10^^ s^^, so that we can identify the fre- 
quency dependent quantity a = cr{uj) with a constant, 
equal to the static conductivity cTstat , only for frequencies 
1/ = uj/2tt less than about 10 GHz. Our frequencies are 
obviously higher: a vacuum wavelength of Ao — 1.06 /im 
corresponds to ujq — 1.78 x 10^^ s"-^. Thus our a{ui) is an 
effective conductivity, defined in accordance with the low- 
frequency relation (2), but no longer identifiable with the 
static conductivity Ustat- The actual value of a{Lo) has to 
be evaluated in each case. 

Take water as an example, for which n — 1.33 + 5.0 x 
10~^ i for the actual wavelength Writing in general 
the complex refractive index as n = n' -\-in" we get, when 
a/eoLo <^ 1, cr/eQLu — 2n'n" , which in our case amounts 
to only 1.33 x 10^^, corresponding to 



(T = 0.21 S/r 



(5) 



This can be compared with the static conductivity for 
pure water, which is much smaller, CTstat = 2 x lO"** S/m 
0. It is notable that the value of ct is real. 

Increasing the incident wavelength somewhat, the ab- 
sorptive effect can be enhanced. Thus, if we consider 
glass, and choose Aq = 8.0 /xm, we have n = 0.41H-0.323 i 
ll^. With w = 2.36 X lO^^^s^i in this case we get for the 
effective conductivity quite an appreciable value: 



a = 550 S/m. 



(6) 



(Actually the parameter u jt^io is as high as 0.265 in this 
case, and can hardly be assumed to be a "small" param- 
eter any longer if Aq is increased further.) 

We assume now that the incoming beam in medium 
1 is a plane wave, such that its electric and magnetic 
field amplitudes are related through = y^co/JiQ Eq. 
We introduce the non-dimensional wave number, a, in 
medium 2, and the relative radial coordinate f: 



a 



— n2U>a/c, f = r/a, 



(7) 



and make use of the standard spherical coordinates 
(r, 9, if) centered at the origin of the sphere. Expansions 
for the incident electric field E^') = {Ei'\ e'^\ eI^^ ) , the 
scattered field = (E:^''', £;^"\ £;^"^), and the inter- 
nal field in the sphere, = {eI'"\e^^\e''^^), are 
for convenience collected in Appendix A together with 
the analogous expansions of the magnetic field. In these 
expressions, we have imposed the electromagnetic bound- 
ary conditions at the spherical surface, r = a. We will 



need use of the Riccati-Bcsscl functions ■)/'; and ^ 



(1). 



(8) 



Here, ji and hi are the spherical Bessel and Hankel func- 
tions, satisfying the Wronskian Wlipi^^i^^} — i. There 
are three groups of expansion coefficients: {Aim, Bim} 
for the incident field, {a;m,^/m} for the scattered field, 
and {cirmdim} for the internal field. As for the first- 
mentioned coefficients, they are determined by inverting 
the formulas IpTTjl and for Ei'^ and H^''^: 



Aim — 
Blrn = 



{b/af 



Eol{l + l)Mk2b)jn 
{b/af 



Ei^\b,e,ip)Y:mie,ip)dn, (9) 



HoKl + l)Mk2b)Jn 



H^\b,e,^)Y*,{e,v)dn,{iQ) 



where d£L — sin 9 d9 dep. Here, the angular integrations 
arc taken over a sphere of arbitrary radius b. This free- 
dom of choice for b gives us a calculational advantage. 
The results for Aim and _B;,„ are only related to the inci- 
dent field, and they have to be independent of b. Calcu- 
lational insensitivity with respect to various input values 
for b thus serves as a check of the evaluations of Aim and 
Bim- In the following, we will suppress the writing of the 
time factors exp(— iwt). 

Our expansion procedure follows the conventions of 
Barton et al. . This procedure was also followed 

in Rcf. In the case where b is chosen equal to a, Eqs. 
© and" Uni) agree with Eqs. (A.23) and (A.24) of Rcf. 

In order to calculate Aim and Bim, it is necessary 
to know the incident field components Er''^ = E",*--^-* and 
i/r'' = Hr^^ in region 2. We introduce the parameters 



niuj 
c 



sm^ 9i 



niLU 



■ sm c/i, 7121 



n2 
ni' 



(11) 
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and take the following expressions for the amplitude ra- 
tios, Til = E, 



(2) /p(l) 



of the substrate into account 



and Ti E^'^'^ /e''1^ at the surface 
2n2i cos Ox 



n^i costal 



(12) 



2 cos 6*1 



(13) 



ii/sm Wi 



Here, ii'ii refers to the field component in the plane of 
incidence [p polarization), and E±_ to the field orthogonal 
to it (s polarization). For the radial field components, we 
find 



= I J_T||i;|'i^^[sin6'i sin 6* cos (ys - i^^sin^ Oi - cos 6*] + Tl£;[^^ sin 6* sin 93 | x exp[-/3(x + /i) + ^72], (14) 



H^'^ ^ |T_ii7|p[-sin6'isin6'cos(p + iysin2 6li -n2^cos6l] + n2iT||ffy^^ sin 61 sin .^j x exp[-/3(a; + /i) + i7z],(15) 

I 



where 



ijf/^fl'^ =^21^11, (16) 



at the surface of the substrate. (Readers interested in 
background material for our developments can e.g. con- 
sult the articles 0, H, O, , as well as the standard 
texts [l|ll[l3-) Upon insertion into Eqs. ® and (jTUl, 
we find that the coefficients Aim and Bim for the case of 
s (_L) and p (||) polarization are related through 



Im 



(17) 
(18) 



hence, it is sufficient to only consider one polarization di- 
rection in the calculations of the expansion coefhcients, 
and in the following, we will choose the case of p polariza- 
tion. By inserting the expressions for the incident field 
into Eqs. © and (|10|l . we obtain the following expres- 
sions for the expansion coefficients: 



ai [I, m) 



sm9iQi{l, m) 



I SIT? 61 — Jl^i Q2{1, ™) 

n2ai{l, to)T|| e"'''* Qsil, m 



where we have defined 



ai{l, m) 
Qi{l,m) 



21 + 1 {l-m,y. 
Att (/ + m)! 



1/2 



27r(-l) 



m—l 



r/2 



sm 



2^1 cos 



(19) 
(20) 

(21) 



(76 cos 0) 



7r/2 



xP["{coa9)Iim\{u)d0, 
ni 

xPf\cosd)Iim\{u)d9, 
where u — pbsinO, I + m is 



smWcosy< * ^ (76COS6') 







COS ] 








i sin J 



(23) 

(76 COS 6*) 
(24) 

and Im{z) = 



even 
odd 

Jm{iz) is a modified Bessel function. (In Ref. j^, 
Eqs. (46) and (49), were missing factors of 2.) Again, 
the calculated values of aJ'^^ and sf^ have to be numer- 
ically independent of the value chosen for the radius b. 

Once Aim and Bi„i are known, the other sets of coef- 
ficients, {a;„i,5;,„} and {Q„i,d/„i}, can readily be calcu- 
lated using Eqs. HA.19|I - HA.22() . In appendix B, we give 
a table of Aim and Bim up to Zmax — 7 for p polarization. 



III. THE SURFACE FORCE 

Let us, as indicated above, write the total radiation 
force F on the sphere as a sum of two contributions: 



Ttsurf 



Tiabs 



(25) 



where F''"''^ is the result of the force density f acting 
in the surface layer, i.e. Eq. (^), and F'^'''' is caused 
by the absorption in the sphere. By integrating across 
the surface layer, the components Ff™^ of F*'"''^ can be 
related to the surface integral of Maxwell's stress tensor 
Sik on the outside of the sphere: 



;isurf 



SikTikdS, 



(26) 



n being the outward normal (cf. for instance, Refs. 
xP™(cos0) [/|,„„i|(w) + /|„j_|_i| (li)] d6, (22) [Til [l5jl. [This expression assumes that the effect from 
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absorption is negligible, in other words that the size 
of the particle is much smaller than the skin depth 
S = .] If a <C Ao, and the susceptibility of the 

particle is small, the formula for the total force F goes 
over to the conventional expression for dilute particles: 



,eok|£;«| 



(27) 



cf., for instance, Refs. [IE 113, Ell- 

It should be stressed that we are calculating the force 



on the sphere in two steps: First, the surface integral 
(26) is evaluated assuming the refractive index to be real. 
That means, the coefficients Aim and Bim below are cal- 
culated on the basis of a real n^. Thereafter, the ab- 
sorptive part F^'^'^is calculated separately, via the force 
density J x B. 

Using the general expressions for the incident and the 
scattered field, Eqs. (|A.1|I - (|A.12|I . we can re- write this 
surface integral in terms of the complex expansion coef- 
ficients Aira, Bira, and bijn- 



^^^^ 



2-^ 



OO / 



l — \ m—~l 



{l + m + 2){l + m + l) 
(2lTT)(2lT3) 



1/2 



l{l + 2) ( 2722aimClt+l,rn+l + '^2'3^'m^r+l,m+l 



{I -m + l){l -m + 2) 



{2l + l){2l + 3) 

X l{l + 2) ( 2n2ai + i^m~ialm + ^2^l + 'i-,m~l^'lm + ^2^l+l,m-iaim + 26/+i,„i_i5;*^ + bl+i^m-lBi^ 

1/2 



1/2 



(28) 



F: 



'surf 



2„2 



n CXD I 



1=1 



{I -m+ + m + 1) 



1/2 



/(/ -I- 2)3 2n\ai+i^ma*i,n + nlai+i^,nA^„^ + n\Ai+i^raa*ir, 



(2;-M)(2Z + 3) 

2bl+i^rab*im + bl+l^mBfj^ + i3i+i,,„6,*„ -|- n2m{2aimb*im + 0-lmBi„^ + Ai„ibi„^) 

I 



(29) 



We shall not dwell much on these expressions, since our 
main interest, as mentioned above, is the vertical force 
F^^^ associated with the absorption. Consider, however, 
the horizontal surface force component F^^"^ . With an 
incident beam power (in vacuum) of P = 150 mW, dis- 
tributed over a circular cross-sectional area of diameter 
10 fim, we find the magnitude of the Poynting vector 
to be {c/2)eoE^ = 19.0 MW/m"^. Taking the radius of 
the sphere a = 1 /zm, we find the denominator in Eq. 

to be eoE^a'^ = 0.13 pN. If the sphere is made 
of glass (rt3 = 1.50) and the surrounding medium 2 is 
water (712 = 1.33), we then have a = 27ra/A2 = 7.9, 
since the wavelength in water for Nd:YAG laser light is 
A2 = 1.06/1.33 = 0.80 ^m. Using Fig. 5 in Ref. 0, we 
find that is positive, as it should be, and is approx- 
imately equal to 0.010 pN in the case of p polarization. 
Since the Reynolds number is very low, the Stokes drag 
formula D = dn^av is applicable, using ^ = 1.0 x lO"'^ Pa 
s as the dynamic viscosity of water. Putting = D, 

we obtain the drift velocity of the sphere to be w = 0.53 /i 
m/s. This is an order of magnitude agreement with the 



observations of Kawata and Sugiura 0; from their Fig. 
4, one infers that w ~ 1 — 2 ^m/s. This agreement is 
satisfactory; there is no need to carry out a complicated 
calculation to find the absorptive correction to the hori- 
zontal radiation force. In the following, we focus atten- 
tion on the vertical absorptive part of the force. 



IV. THE VERTICAL ABSORPTIVE FORCE 

To calculate the absorptive force, we first replace Eq. 
by the general expression 



f = /9E + JxB - -E^\/e, 



(30) 



where p and J are the external charge and current den- 
sities in the medium. A general derivation of this force 
expression can be found, for instance, in Stratton's book 
[3. Also, in the review article of one of the present au- 
thors E3 , it was found that the force expression is able to 
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describe the outcome of known experiments. We assume 
no external charges to be present, so that p = 0. Writing 
J = (tE where a is the frequency-dependent conductiv- 
ity, assumed to be real (cf. Eq. (5)), we get in complex 
representation the following expression for the absorptive 
force: 

pabs = Ij^y J X WdV = y E X H* dV, (31) 

where the integral is evaluated over the volume of the 
sphere. 

In cases of current interest, we can simplify our cal- 
culation: In view of the assumed smallness of a/eow, we 
replace the fields inside the sphere with those calculated 
when assuming the material to be a perfect dielectric. 
That is, we can replace the complex refractive index 713 
by its real part n^. This facilitates the calculation of the 
volume integrals. The expressions for the internal fields 
are given by Eqs. ljA.13|l - (|A.18(l . in spherical coordinates. 
For notational convenience, we will omit the superscripts 

With f — r I a we can write the vertical absorptive force 

as 



F: 



abs 



07x0^3 
2 



/'2'7r f'T: /•! 

c?^ \ dip smOdO I f^df 
Jo Jq Jo 



(E X H*)^ 



X sin6'cos<y9 +(E x H*)g cos^cosiys — (E x H*)^ siniy9 

(32) 



where we have defined 



p27r pTT pi 

1=1 cos(pdip / sin^ 9d9 / f^dfEgH*, (33) 
Jo Jo Jo 



h 
h 
h 
h 
h 



2,r ,-K ^ 

cos ip dip / sin 
Jo 



f'dfE^Hl (34) 



/ cos ip dp / cos^sinfldfl / r^drE^H* , (35) 
Jo Jo Jo 

- I COS p dip I cos 6* sin Me / p'^dfErH^, (36) 
Jo Jo Jo 

fill fTT /•! 

/ smp)dp} / smOde I f^dfErH*g, (37) 
Jo Jo Jo 

(38) 





sin ip dp I sin 
"'0 



r^drEoH* 



It is now convenient to introduce the operators 



/=i j=i 



du, 



(39) 



m— — l k——j 

We write the two first of the liS as a sum of four terms: 



Ii = iiTaniceoEl'^Iik, i = 1,2. (41) 
fe=i 

and remaining li 's, which are simpler since they are each 
composed of only two terms: 



li ^ inaniceoEQ^^Iik, i = 3,4, 5,6. (42) 

k=l 

where (in the following, I, m, j and k denote summation 
indices as defined in Eqs. H39|l and (001)) 



hi = -nz2C 



'13 



'14 



/21 



'22 



4^[{u)4,'^ {u)M\^^ I kci„J*kChnC,kni {I, m, J, k) 



-7132^2/1 



n2 



ipi {u)'ipj {u)m'iP I mkdimdjkCimCjk'R'sil, m, j, k) j 
Mu)i^,{u)M\-^Hmd 

imC|fcCimCjfc^l(jj I, m) > 



71,32/1 



i^'Mi;'^ iu)M[^'^ \ mci„,d*kChnCjkni (j, k, I, m) 



i;'i{u)i;j{u)M\^U mkci^c*^.Ci.,nCjkT^3{l,m,j,k) 



(43) 
(44) 
(45) 
(46) 
(47) 
(48) 
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^23 — 



'24 



^31 



'32 



141 



142 



hi 



152 



"2 



"32 

1 

n2 

"32 



,m,j, k) 



,m,j,k) 



+ l)M^i^H ■mcimd*kCimCjk'R-4{l, m,j, k) 



C 



u 



,m,j, k) 



^^^'jiu)l{l + 1)M\P\ kci„4*kClraCjkTl4{l,m,j, k) 



2 r 



'"32'C 



^j{u)l{l + 1)M\^^ l^ir,tC*i.ChnCjk'Tlb{j, k,l,m) 

^^^r{u)l{l + l)M'i;^ \ Cl^d*,,ClrnC,Mj, k, I, m) 



n32n2>C 



+ l)M'i;Ukcir,,C*^ClmCjkn3{l,m,J, k) 



hi = n-iiC. 



h2 — 



"2 



^^'^"V j(^')jO' + \ mdimd*kCimCjkTi3{l,m,j,k) 



(49) 
(50) 
(51) 
(52) 
(53) 
(54) 
(55) 
(56) 
(57) 
(58) 



r 



with 



de 



^ „ , ' ■ 2. dPr [cos 9) dPNcos 



d9 



Pf{coB6)de, (59) 
de, (60) 



de 



n3{l,m,j,k) = / P/" (cos 61) PNcose) de, 



(61) 



7^4(^, TO, j, fc) = / cos e Pi"' (cos 6i) P/ (cos 61) d6i, (62) 

dP"^( cose) 

n^{l,rn,j,k) = / cos 61 sin 6* — -Pf{cose)de, (63) 



and the definition 



'(2/ + 1)(Z-to)! 



(64) 



47r(? + TO)! 

Altogether, writing the non-dimensional vertical absorp- 
tive force as Q^^, where 

^abs 



Q 



abs 



we get 



ni2a 



4 



2 4 



6 2 



1=1 fe=l 



i=3 fe=l 



(65) 



(66) 



Note that the prcfactor in the above equations, as well as 
the terms hk, are non-dimensional. The expression H66|) 
is our main result. 



V. RESULTS 
A. Numerics 

We implemented an adaptive Gauss-Kronrod rule to 
calculate the integrals in Eqs. (HU - and (gSJ - 
(|58|l . By employing the symmetry properties of the inte- 
grals over the Legendre polynomials (Eqs. lf3^ - ^^), 
we were able to greatly reduce the number of integrals to 
be computed. To make sure that the sum over I and j in 
Eq. H39|l converged properly for our chosen cut-ofT value, 
we plotted the total absorptive force on the sphere versus 
increasing cut-off, /max, in I (we used the same cut-off for 
I and j). In Fig. 12 we give an example for p polarization 
using several values of a. The total force is normalized 
to the Zniax — 30 result; all other parameters are given 
in the figure caption. This figure clearly suggests that, 
for a given value of a, there is a narrow range in I which 
accounts for most of the absorptive force. Moreover, it 
is seen that this range of important /-values moves to 
higher / and broadens with increasing a. In the inset, we 
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FIG. 2: Plot of the total absorptive force versus increasing 
cut-off in Eq. ^ for 4 values of a: 2.5, 5.0, 7.5, and 10.0, 
with ni = 1.75, n2 = 1.33, and ns = 1.50. The inset shows 
the absolute value of the difference in the total force using I 
and (/ — 1) respectively, as cut-offs. 



plot on a linear-log scale the absolute value of the dif- 
ference in the total force using I and — 1) respectively, 
as cut-ofFs. The inset indicates that contributions to the 
absorptive force from large ^-values (large with respect 
to the Z-value of the peak) are exponentially suppressed. 
Hence, for the range of values for a, ni, n2, and na stud- 
ied in this article, we need only consider contributions to 
the force from I < 30. 

It is of interest to point out that the main features 
of Fig. 121 can be understood from a physical point of 
view. (This kind of argument is frequently made use of 
in connection with Casimir-related calculations; cf. for 
instance, Ref. 19] ). The most significant angular mo- 
menta are those that are of the same order of magnitude 
as, or are somewhat smaller than, the angular momen- 
tum corresponding to grazing incidence on the sphere. 
Consider now a photon with energy huj and momentum 
Hk2 that just touches the surface of the sphere. Its angu- 
lar momentum is hk2a. Setting the angular momentum 
equal to hi, we thus obtain I — — a. Thus, the most 
significant values of I should be expected to be / < a. 
This is seen to be in accordance with the calculated re- 
sults in Fig. 121 In view of the crudeness of the argument, 
the agreement is actually better than we had reason to 
expect. 



B. Discussion of the figures 

Illustrative examples of the behavior of the non- 
dimensional vertical absorptive force are shown in Figs. 

The angle of incidence is 9i — 51° everywhere. 
In the case oi p polarization, Eq = whereas for s 



FIG. 3: Non-dimensional vertical absorptive force = 
Fx^^ /{^oEaa^) versus particle size parameter a = k2a for the 
set of refractive indices shown. The two states of polarization 
for the incident beam are distinguished. Here, as well as in 
the subsequent figures, 6i — 51°. Also, all figures refer to 
sphere resting on the plate, i.e. to h = a. The figure is scaled 
against the conductivity a of the sphere. 



polarization, Eq — E^^\ When the incident field is taken 
to be a plane wave, the amplitudes and e''^^ are 
constants. It is to be kept in mind that, in order to lift 
the sphere from the plate, the positive absorptive force 
has to overcome not only gravity, but also the negative 
surface force FJ"'^^ that was calculated in Ref. Q for the 
same values of {nr, n2, ns} as we consider here. The cases 
treated in the two papers are thus directly comparable. 

Similarly to Ref. 2], all curves in Figs. 0- [3 refer to 
the case of contact between sphere and plate, i.e. to 
h — a. From Eqs. (|19() and (|20|l . it is seen that Aim and 
Bim, contain exp(— /3/i) as a common factor. The same 
property is carried over to the coefficients q„i and dim, 
according to Eqs. I|A.21I) and IIA.221) . As the absorptive 
force is quadratic in the last-mentioned coefficients, it 
follows that Q^'^ at an arbitrary height h is equal to 
exp[—2(3{h — a)] times the value that can be read off from 
Figs.|31-[S1 No figures need to be worked out to show 
the dependence of the force with height above the plate. 

In order to interpret the figures physically, it is helpful 
to consider a concrete example. Let us relate our discus- 
sion essentially to the situation considered in Section 3: 
incident plane wave of power P = 150 mW, correspond- 
ing to eQEQO^ — 0.13 pN when the radius of the sphere is 
a = 1 nm. If the density of the sphere is 2.4 g/cm"^ (the 
same density as for glass), the weight of it is mg — 0.10 
pN. Moreover, 713=!. 50, n2— 1.33 (water), whereas the 
size parameter is a = 7.9 for Nd:YAG laser light. Let 
us look for the necessary value of a in order to satisfy 
the condition for elevation of the sphere, just above the 
plate. The condition is F^^^ = mg + \F^™^. Expressed 
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FIG. 4: Same as for Fig. but with a higher (dominant real 

part of) refractive index (713 = 1.60) in the sphere. FIG. 5: Same as for Fig. but with vacuum (na = 1) sur- 

rounding the sphere. 



non-dimensionally, 

Qf^ = ^f^ + iQr'i- (67) 

From Fig. |31 we read off, when choosing for definiteness 
the case of p polarization, Q'^^ /{lO^'^a) — 1.7, whereas 
from Fig. 4 in we read off Q™"^^ = —0.35 for the same 
value of a. This leads to a relatively high value of cr, 
about 6600 S/m (S= fi^^). In our example, the influence 
of the weight of the sphere is relatively large because the 
incident power is so moderate. If we increase P to 1 watt, 
we get eo-Egfl^ = 0.87 pN, resulting in a somewhat lower 
value, a « 2700 S/m. 

For glass, at optical wavelengths, the effective con- 
ductivity a is very small. (The static conductivity is 
extremely small, Cstat ^ 10""'^^ S/m.) Even for higher 
wavelengths such as Xq = 8.0 /^m, the value a — 550 S/m 
given in Eq. (6) is insufficient to account for the lift- 
ing force. For water, at optical wavelengths, as we have 
seen in Eq. (5), the effective conductivity is even smaller. 
We can thus immediately conclude that absorption is not 
the physical reason for the elevation of the sphere in the 
Kawata-Sugiura experiment. It is likely that there were 
surfactants on the spheres in this experiment, and hence, 
that a theoretical description of the kind given in Ref. 
PI is most appropriate in this case. For other materials 
where the conductivities are higher, for instance for car- 
bon (cstat — 0.77 X 10'"' S/m), the absorption-generated 
elevation of microspheres should be quite possible physi- 
cally. One must however be aware of the algebraic restric- 
tion made in the present paper to simplify calculations: 
our calculation relied on the property that a/eoto <^ 1 
(cf. Eq. (0)). This restriction made it possible to take 
the refractive index 713 to be real, inside the integrals in 
Section 4. Since cqcv « 10^ S/m, we can no longer expect 
high accuracy from the formalism when cr 3> 10^ S/m, 



at optical frequencies. 

Consider finally the remaining figures: Fig. ^ shows 
that when the refractive index in the sphere increases 
from 1.50 to 1.60, the absorptive force increases. This is 
as we would expect. Moreover, if the difference between 
interior and exterior indices increases, the Riccati-Bessel 
functions cause the force to become more oscillatory in 
character. The changes are clearly shown, if we com- 
pare Fig. 131 with Fig. El Increasing differences between 
the refractive indices effectively mean enhancing the the 
geometrical-optics properties of the system. There is also 
a dependence of the force upon the two different polariza- 
tion states of the beam. No simple physical explanation 
of this dependence seems to exist. 



VI. CONCLUSIONS 

We can summarize as follows: 

1. The strategy of the above calculation has been to 
identify the evanescent field with the incident field on 
the sphere. The theory makes use of wave optics, and is 
thus applicable even when the wavelength is of the same 
magnitude as the sphere diameter. 

2. The main novel development of the present paper 
in comparison with Ref. '5| is the calculation in Section 
4 of the vertical absorptive force, assuming the refractive 
index 77.3 in the sphere to be complex. The calculation 
rests upon the approximation that a/eauj ^ 1, so that 
the fields inside the sphere can be evaluated as if the 
material were a perfect dielectric. The surface force is 
thus calculated from Eq. (26) assuming = to be 
real. The absorptive force is calculated from Eq. (30). 

3. The results of Figs. - [31 show that in order to 
account for elevation of the micrometer-sphere, the con- 
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ductivity cr must be at least as large as about 10^ S/m. 
This applies, for instance, to media like carbon. 

4. For latex or glass spheres, used in the Kawata- 
Sugiura experiment 1], the absorptive force seems to be 
too weak to account for the elevation force. Most prob- 
ably, the observed elevation force was due largely to im- 
purities or a film on the surface of the sphere, making it 
partly conducting. Cf. also Ref. The absorption is 
sensitive with respect to impurities, and it is at present 
difficult to estimate to which degree the absorption is 
important in practice. 

5. As a final general remark: the subject of radiation 
pressure on particles in the evanescent field is a many- 
facetted phenomenon which is not yet fully understood, 
and it is a very active subject. In addition to the ref- 
erences given above, the reader is referred also to recent 
papers in which the dipole method is made use of, in 
order to calculate the force of an evanescent wave over 
dielectric and metallic surfaces ^20, 21]. 
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APPENDIX A: EXPANSIONS OF THE 
ELECTROMAGNETIC FIELDS 



For convenience, we summarize the expansions for the 
incident, (i), the scattered, (s), and the internal, (w), 
field components. As mentioned in the main text, sub- 
script 1 refers to the non-absorbing substrate, subscript 
2 to the non-absorbing medium around the sphere (re- 
fractive index ^2), and subscript 3 to the absorbing 
sphere whose complex refractive index is n^. The rela- 
tive refractive index is = ^3/''^2- Furthermore, Eq 
and Hq refer to the incoming plane wave in the sub- 
strate, the non-dimensional wave number and distance 
is a = = n2U>a/c, and f = r/a respectively. 

The formulas below are general. As mentioned in the 
main text, in the practical calculations of the coefficients 
Aim and Bim we took however the refractive index 713 to 
be real. 
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EE 
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H. 



Hi-) 
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, -x^^^m , ^(1). 



' sm fc* '■ do 
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Internal field 



OO / 



/^l m— — l 

j-i GO / 

l — l m— — l 

4'^'=^E E 



«32 Qm Wi[n32ar)—- dim Wi{m2ar)- 



86 n2 



sin 9 



/ — I m— — / 

00 / 



imn32Cim-4''i{n32af) - —dimipi (n32af) 

sm (7 n2 uu 



m— — l 
. . 00 I 



/^l m—~l 
00 / 



at^ sm ( 



Z — Im— — / 



imn32dhnip'i {n32af) ■^^ + in2nl2Cim'ipi {n32af) 

sm t' at' 



(A.13) 
(A.14) 
(A.15) 
(A.16) 
(A.17) 
(A.18) 



Everywhere, primes denote differentiation with respect 
to the whole argument. 

Relations between coefficients 



i'i{n32a)£,l^\a)-n32tp'i{n32a)^l^\a) 
Here the Wronslcian is Wlipij^l^^} = i. 



Blrn- (A.22) 



i>'i{n32a)-ipi{a) - n32'ipi{n32a)^p[{a) 
n32il'i{n32a)£,l^\a) ~ip[{n32a)^l^\a) 



A. 



Im 1 



, _ n32ip'i{n32a)ipi{a) - ^i{n32a)ip'i{a) „ 



V'i("'32a)CP^(")-"32'0;'("-32a)CP\a) 
i 

nl2-ij)i{n32a)s\^\a)-n32-ip[{n32a)s!i^\a 



(A.19) 
(A.20) 



■Aim, (A.21) 



APPENDIX B: TABLE OF Aim AND B;™ 

We include a table of the numerical values for the real 
and imaginary parts of Aim and Bim for an evanescent 
field. The coefficients are calculated for p polarization, 
and we have set a = 1 and e~^^^ = 1 in Eqs. l(T^ and if^ . 
For space considerations, we only give the coefficients up 
to / = 7. 
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